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Abstract. We prove estimates on the Lebesgue constants of the hyperinter- 
polation operator for functions on the unit ball B d C R d , with respect to 
Gegenbauer weight functions, (1 — |jx|| 2 )'' -1 '' 2 . The relationship between or- 
thogonal polynomials on the sphere and ball is exploited to achieve this result, 
which provides an improvement on known estimates of the Lebesgue constant 
for hyperinterpolation operators on B 2 . 



1. Introduction 

Hyperinterpolation was introduced by Sloan in [7] as an approximation technique 
which preserves a classical inequality of Erdos and Turan, concerning interpolation 
polynomials in a univariate setting, to higher dimensions. Combining two common 
approximation techniques, hyperinterpolation uses the Fourier orthogonal expan- 
sion of a function in orthogonal polynomials, but approximates the integrals used in 
the expansion by means of a cubature formula with positive weights. Hence, hyper- 
interpolation provides a polynomial approximation which relies only on a discrete 
set of data. 

Much attention has been given to hyperinterpolation on spheres and balls. On 
spheres, Le Gia and Sloan [4] studied the Lebesgue constant of the hyperinterpola- 
tion operator, showing that Lebesgue constant grew on the order of the Lebesgue 
constant of the minimal projection operator on C(S d ~ 1 ), subject to certain "reg- 
ularity" conditions on the cubature formulae. Reimer was able to independently 
prove in [6] the same estimate obtained by Le Gia and Sloan, and further, that 
the "regularity" condition imposed on the cubature formulae was not needed, and 
that the result held for any positive cubature formula. Specifically, if C n is a hy- 
perinterpolation operator on the space of continuous functions on S d ^ 1 , which is 
a projection onto the space of polynomials of degree n, Le Gia and Sloan, and 
Reimer, were able to prove there exist postive constants bd and Cd satisfying 

W d - 2)/2 <||/: n ||< c ^- 2 )/ 2 . 

Similar estimates were obtained for Sobolev spaces on the sphere in [3|. 

Regarding hyperinterpolation on the ball, Atkinson, Chien, and Hansen pQ in- 
vestigated hyperinterpolation on continuous functions on the unit disk B 2 in R 2 . 
An esimate for the Lebesgue constant for a hyperinterpolation operator using a 
particular cubature formula was found. The cubature formula was obtained by 
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combining quadrature formulas on a closed interval and a circle, and the estimate 
on the Lebesgue constant obtained for the hyperinterpolation operator was 

H4.ll < cnlog(n+l) 

for some positive constant c. 

In this paper, we improve the estimate obtained in [TJ, and extend the results to 
general hyperinterpolation formulas stemming from cubature formulas with positive 
weights on balls in higher dimensions, where the cubature formulas may also be with 
respect to certain Gegenbauer weight functions. Our results rely on the results on 
the sphere obtained in [6], and the connection between orthogonal polynomials on 
the sphere and the ball, which were developed by Xu. 

In Section 2 of this paper, we will cover the background information needed for 
the proof of our theorem. In Section 3, we will state and prove our main theorem 
on the Lebesgue constants of the hyperinterpolation operators. 

2. Orthogonal Polynomials and Cubature Formulae 

2.1. Orthogonal Polynomials on B d and S . The structural interplay between 
orthogonal polynomials on the unit ball B d and S d are essential in our result, so we 
review these results here. It is assumed that the reader has a basic understanding 
of orthogonal polynomials; the standard references are [9] for univariate and [2] for 
multivariate orthogonal polynomials. 

We denote the space of polynomials in d variables with total degree less than or 
equal to n by LT^. We define V d to be the space of polynomials of total degree n in 
d variables. 

We denote the space of orthogonal polynomials of total degree n in d variables 
with respect to the weight function to on a set f2 by V d (uj). We will be particularly 
interested in orthogonal polynomials on B d , with respect to the Gegenbauer weight 
function, w M (x) = (1 - 1 1 1 2 ) jLi_ 1/2 for /! > -1/2. 

The dimension of the space V^w^) is r d — Given a orthonormal basis 

of V^(w M ), {P£ v P£ 2 , . . . , P^ r i}, the reproducing kernel of V d , P„(w M ; x, y), is 
the function defined by 

(2-1) Pn(w M ;x,y) = £ P n ^(x)P£.(y). 

It earns its name because it "reproduces" orthogonal polynomials of degree n: 

/( x ) = / Pn(w M ;x,y)/(y)w Ai (y)dy 
J B d 

for all / £ V d . The reproducing kernel of the space II„ with respect to w M is the 
function K n (w Jtl ;x, y), defined by 

n 

K„(w M ;x,y) = ^ Pj(w^; x, y) 

3=0 

and satisfying 

(2.2) /(*)=/ K, i (w M ;x,y)/(y)w M (y)dy 

J B d 

for all / e n£. 
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We now turn to orthogonal polynomials on the sphere. Orthogonal polyno- 
mial structure on the unit sphere differs somewhat from the orthogonal polynomial 
structures observed in Euclidean regions with non-zero measure, due to the fact 
that unit sphere is a algebraic surface. The orthogonal polynomials on the unit 
sphere with respect to the surface measure dm (du>d for S d when dimension requires 
specification) are known as spherical harmonics, which are harmonic, homogeneous 
polynomials restricted to the unit sphere. It turns out that harmonic, homogenuous 
polynomials of different degrees in d + 1 variables are orthogonal with respect to 
integration over S d . We denote the space of homogeneous, harmonic polynomials of 
degree n in d+1 variables by When restricted to the unit sphere, the space of 

polynomials of degree n can be decomposed orthogonally into the subspaces Hf. , 
where we are considering the polynomials to be restricted to the unit sphere. 

Spherical harmonics can be generalized in order to include the possibility of 
weight functions; for example, a study of spherical harmonics with respect to weight 
functions invariant under reflection groups is given in Chapter 5 of [2]. When 
considering spherical harmonics with respect to a given weight function H , we use 
the notation U d+1 (H). 

A reproducing kernel for spherical harmonics also exists. We may take an or- 
thonormal basis for the space of spherical harmonics of degree n to obtain the 
reproducing kernel for that space, in the same way as (|2.1j) . We denote the repro- 
ducing kernel of the space of spherical harmonics of degree n with respect to the 
weight function H by P„(i?;x, y), and it satisfies 



where ojd is the standard surface measure for S , whenever / is a spherical har- 
monic of degree n. We may also obtain a reproducing kernel for the space of all 
polynomials, restricted to the unit sphere, of degree less than or equal to n, K n , by 
letting K n (H;x,y) = H£ =0 -Pfc(-H";x,y); K n then satisfies 



whenever / is a polynomial of degree less than or equal to n and x G S . 

2.2. Relationship between Orthogonal Polynomials on Spheres and Balls. 

The relationships between orthogonal polynomials on spheres and balls is very 
rich. We will focus on the relationship between reproducing kernels on spheres and 
balls. This subsection outlines this relationship, which is given detailed treatment 
in Section 3.8 in [2]. We first provide some definitions that are needed to describe 
the relationship we will exploit. 

Definition 2.1. A function H is centrally symmetric if -ff(x) = H(—x) for all x 
in its domain. 

Definition 2.2. A function H is positively homogeneous of order a if, for all x in 
its domain and all t > 0, H(tx) = t a H(x). 

Definition 2.3. A function H is admissible on R d + m + 1 if H(x) = i/i(xi)i/ 2 (x 2 ), 
where xi £ R d , x 2 G IR™" 1 " 1 , Hi is centrally symmetric, and H2 is positively homo- 
geneous. 



(2.3) 




/(x)=/ K n (H;x,y)f(y)H(y)du Jd (y) 



Js d 
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With these definitions, we are able to state the following theorem, which will be 
crucial in our results. 

Theorem 2.4 (Thm 3.8.12, [2]). Let d and m be positive integers. Suppose that 
H is an admissible weight function on as in Definition \2.3l where Hi is 

positively homogeneous of order It. Let W^(y) = -Hi(y)w T+m / 2 (y) &e defined on 
B d . Then the reproducing kernel of V n on B d with respect to W^ 1 , P„(W^;-, •), 
is related to the reproducing kernel of spherical harmonics on S d+m with respect to 
the weight function H , P n (H; •, •) by the formula 

(2.4) P„OT;x,y)= f P n (H; (x,x c ), (y, y/l - \\y\\ 2 n))H 2 ( V )d^ m (n), 

Js m 

where x c £ R m+1 is any point satisfying (x, x c ) G S d+m+1 . In particular, if Hi = 
Hi = 1, then 

(2.5) P n (w m/2 ;x,y)= / P n (l; (x, x c ), (y, y/l - \\y\\^))duj m ( V ). 

We may sum both sides of the formula in (12. A} to obtain the following corollary. 
Corollary 2.5. Under the assumptions of Theorem \2.4\ we have 

(2.6) K„TO;x,y)=/ K n (H; (x,x c ), (y, y/l - \\y\\ 3 v))H 2 (T,)dLJ m (v), 
where x c has the same meaning as in Theorem \2.4\ If Hi = Hi = 1, then 

(2.7) K„(w m/2 ;x,y) = f K n (l; (x,x c ), (y, y/l - \\y\\*r,))d«>m(v)- 

2.3. Cubature and Hyperinterpolation. Using the reproducing kernels, one 
is able to obtain the best polynomial approximation in the L 2 norm to a given 
function. Hyperinterpolation uses a cubature formula to approximate the integrals 
required to obtain this approximation. 

Definition 2.6. Given a subset fl o/M d with weight function u>, a cubature formula 
of degree n is a set of positive weights (A Q ) C M + and nodes (x a ) C where a 
belongs to some index set A, satisfying 

^2 ^af{x a ) = / f(x)uj(x)dx 

whenever f is a polynomial of degree less than or egual to n. 

We will make use of an arbitrary cubature on S d in our proof. Many examples 
of such cubature formulae are known. For examples, see [S], [5] and Theorem 5.3 
in [TO]. 

Given a cubature formula on B d or S d , the hyperinterpolation operator is ob- 
tained by approximating the integral in the formula (I2.2[) or (|2.3[) with an appro- 
priate cubature formula. Specifically, taking a cubature formula of degree n, then 
the linear operator 

£„(/)(x) := ^2 X a K n (x a ,x)f(x a ) 
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satisfies 

£„(/)(*) = /(x) 

for all / G Il„. The operator C n is the hyperinterpolation operator of degree n. 

3. Main Results 

The main result of this paper is the following theorem, concerning the Lebesgue 
constant of C n , which is defined by 

||£„|| = sup{||£„/|| c(B , ) : / G C(B d ), \\f\\ C (B«) < !}• 

Theorem 3.1. Let Lm be the hyperinterpolation operator arising from a cubature 
formula with positive weights on the ball B d , d > 1, with respect to the Gegenbauer 
weight function w M (x) = (1 — 1 1 >c 1 1 2 ) ^ !/ 2 ^ where pi is a positive half-integer. Then 
there exist positive constants b d ^ and c d ^ satisfying 

b d .M d - 1),2+ "<\\^ n \\<c d .^ d - l )l^. 



Integral in the proof of our theorem is the following result on hyperinterpolation 
on spheres proven by Reimer in [6] . 

Theorem 3.2 (Theorem 1, [6]). If C^n is the hyperinterpolation operator arising 
from a cubature formula of degree 2n with positive weights on 5 ld_1 , then the norm 
of C,2n as a projections operator from C(S d ~ 1 ) to LT 2 „ satisfies 

a d n^ 2 <||A„||<W d - 2)/2 , 

for positive constants a d , b d . 

Before we prove Theorem 13- 11 we first provide two lemmas that illustrate the 
relationship betwee cubature formulae on B d with respect to the weight w m / 2 and 
cubature formulae on S d+m . Both of these lemmas rely on the integration formula 
below. 

Lemma 3.3 (Lemma 3.8.9, [2]). Suppose that d and m are positive integers. Then 



(3.8) / f(y)du(y)= w m/2 (x) / /(x lV /l- ||x|| 2 77)<M7,)dx. 

Jgd+m J B d J Sm 

Our first lemma shows that we may use a cubature formula on S d+m to obtain 
a cubature formula on B d with respect to w m / 2 - 

Lemma 3.4. Suppose that ^2 a< zA^af{x a ) is a cubature formula of degree n on 
S d+m . This cubature formula is also of degreen n on B d with respect to the weight 
function w m / 2 , in the sense that 



/ /(x) w m/2 dx = y2 \ a f(x a ) 



whenever f is a polynomial of degree less than or equal to n, where f is a function 
defined on S d+m , but only depending on the first d variables, and satisfying 

where \ S m \ is the surface measure of S m . 
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Proof. First observe that 

/ d ^W Wm / 2 W dx = y y d /(x,i7\/l -x 2 )w m/2 (x) dxdu} m (v) 



1 



15" 



/(y)rfw d+m (y), 

where we have used Lemma 13.31 for the second equality. Applying the cubaturc 
formula to the last integral immediately gives the result. □ 

The next lemma illustrates how a cubature formula on B d with respect to w m / 2 
can be used to obtain a cubature formula on S d+m . 

Lemma 3.5. Suppose that ^2 aeA X Q f(x a ) is a cubature formula of degree n on B d 
with respect to w m / 2 (x y ). Let X^es v &f{tf)) ^ e a cubature formula of degree n on 
S m . Then the cubature formula 

aGA,0eB 

is a cubature formula of degree n on S d+m . 

Proof. Let / be a polynomial of degree n on M. d+m+1 . Note that, in (|3.8p . any 
odd powers of y/l — ||x|| 2 will vanish after integrating over S m , since an odd power 
of y/1 — ||x|| 2 in the integrand must be multiplied by at least one odd power of a 
component of 77, and any odd component of n will result in zero upon integrating. 
Hence, after integrating over S m , we are left with a polynomial of degree less than 
or equal to n in x. 
It follows that 

^2 ^2 ^/3f(x a , y/l- \\x a \\Hp) = K f{x a , V 1 - \\x a \\ 2 r))duj m (ri) 



/(x,Vl-||x|| 2 ^ m (?7) 

B d JS m 

x w„ l/2 (x)dx 

/(y)wd+m(y), 

Sd+m 

so that the extension is a degree n cubature on S d+m . □ 

We now prove Theorem 13.11 

Proof. We first use the standard inequality for estimating Lebesgue constants, and 
then apply Corollary 12.51 to obtain 



||^2n|| < sup X a \K n (w M ;x,x a )| 
^ Bd a eA 

< sup ^ X a I K n ( 1; (x,x c ), {x a ,\Jl - \\x a \\ 2 n)\ du} m (rj) , 
^ Bd a eA J S m V ' 

where we have let /1 = 2m. Let vpfitp) be any cubature formula of degree 
2n for S m . Since the integrand is a polynomial of degree n in 77, we may apply a 
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cubature formula of degree 2n to the integral over S m and obtain 



|£ 2 n|| < SUp V A c 
xGS d 



), (2 a , x/l - lk«H 2 ^) 



(3.9) 



< sup y~] V" \ a V0 K n ( l;y, (x Q , \/i - Hxall 2 ^) J 

v pM+m ^ " ^ ' \ / 



By Lemma 1331 the cubature formula J2 a 8 ^a^Bf(x a , yX — II^qII 2 ^) is a cubature 
formula of degree 2n on S d+m . In the proof of Theorem [321 it was proven that the 
expression in (|3.9|) is less than &d+ m (2n)( d+m-1 )/ 2 ; hence we immediately obtain 

\\C 2n \\ < c d+m nC d+m - 1 )/ 2 

for a positive constant Cd+ m . 

To finish the proof, we again let vpf{tp) be a cubature formula of degree 2n 
on S m , and observe that if / 6 C(B d ) with H/llcfS 1 *) < 1) then we may define / 
on S d+m as in Lemma 13^41 and obtain 



^ A Q /(x Q ) = ^ \ a Vfjf(x a , y/l- \\x a \\Hp). 



a£A a, 8 



However, by Lemma 13.51 the cubature formula on the right of the equality is a 
cubature formula of degree 2n on S d+m , and / is in C(S d+m ) and has uniform 
norm less than or equal to |5 m | _1 . Hence, by Theorem 13.21 there is a positive 
constant bd, m satisying 



(d+m-l)/2 



> b d , m n 

which completes the proof of Theorem 13.11 □ 

An immediate result of this theorem is an improvement on the Lebesgue constant 
estimate in [1]. 

Corollary 3.6. Let C n be a hyperinterpolation operator of degree n on B 2 with 
respect to the standard Lebesgue measure. Then there exist positive constants b and 
c satisfying bn < \\C n \\ < cn. 
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